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We extend the analysis of photon coincidence spectroscopy beyond bichromatic excitation and 
two-photon coincidence detection to include multichromatic excitation and multiphoton coincidence 
detection. Trichromatic excitation and three-photon coincidence spectroscopy are studied in detail, 
and we identify an observable signature of a triple resonance in an atom-cavity system. 
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I. INTRODUCTION 

Cavity quantum electrodynamics (CQED) JlJ in the 
optical domain is rapidly progressing: advances in atom 
cooling methods, as well as improved optical cavities 
which allow large dipole coupling strengths, are leading 
experiments into new frontiers of research. Single-atom 
experiments are now possible , and trapping of atoms 
in optical cavities should soon be feasible Q. Quantum 
effects and quantitative testing of theoretical models for 
the CQED system can be performed better than ever and 
perhaps directed to certain applications such as quan- 
tum logic gates H. Exciting developments are also tak- 
ing place in the microwave domain |^-^| , but here we are 
concerned with photon coincidence measurements, which 
are performed only in the optical domain. 

The method of photon coincidence spectroscopy (PCS) 
has been introduced as a means to study the spectrum of 
the combined atom-cavity system p|-[T^1 , but this method 
was restricted to probing only the first couplet of the 
nonlinear regime of the Jaynes-Cummings (JC) spec- 
trum (see Fig.[l]). Here we generalize the method of 
photon coincidence spectroscopy to show how probing of 
higher levels of the spectrum can be performed, and we 
show that three-photon coincidence spectroscopy (3PCS) 
could yield a signature of the second couplet in the non- 
linear regime of the JC spectrum ]l3fl , thus enabling di- 
rect, unambiguous probing of the quantum features of a 
single atom in an optical cavity. 

Photon coincidence spectroscopy is necessary to probe 
quantum features of the atomic CQED system because 
practical difficulties limit the efficacy of other techniques. 
The major difficulties in the optical regime include the 
width of the atomic beam traversing the cavity, motion of 
the atom through the cavity, and the interruption of the 
Rabi oscillation numerous times during passage through 
the cavity as well as a fluctuating atomic number. How- 
ever, for sufficiently slow moving atoms Jlpf , the atoms 
can be regarded as being essentially stationary, and the 
motion and spread of the atoms are responsible, then, for 
an inhomogeneous broadening of the spectral peaks. Fur- 
thermore, at low densities, single-atom effects dominate 
over multi-atom process, and the JC model provides an 
excellent description. Photon coincidence spectroscopy 
was then devised as a way to probe the interesting quan- 



tum features in the presence of unavoidable inhomoge- 
neous broadening. 

Until now, only two-photon coincidence spectroscopy 
(2PCS) has been studied, both as quantum trajectory 
simulations ||Io| and analytic, continued-fraction methods 
]TT] , p~2[ . However, 2PCS is only useful for probing the 
second couplet of the JC ladder, or, equivalcntly, the first 
couplet of the nonlinear regime (see Fig.|I|). Our aim 
here is to propose multi-photon, or TV-photon coincidence 
spectroscopy (NPCS) as a method for probing higher- 
level states and to show explicitly how 3PCS would work 
and its feasibility. 

II. DYNAMICS 

The JC Hamiltonian for the atom-cavity system 

H{g) = fiwAf + ihgA, (2.1) 

where w is both the atomic transition frequency and cav- 
ity resonance frequency, g is the dipole coupling strength, 

Af = <7 3 + fit a + 1/2 (2.2) 

is the "excitation number" operator, and 

A = a) a- - aa + = -A\ (2.3) 

describes the isolated atom-cavity system. The coupling 
strength g is not actually a fixed constant in the calcu- 
lations: instead a distribution P(g) is constructed which 
accounts for the variability of atomic position within the 
cavity mode and results in inhomogeneous broadening 
JuJITlJ. The master equation 

p= [H(g),p]/ih+ [£{t)a+-£*{t)a-,p] 
+ (7j/2)(2<7_ y o<7 + - o+o-p - pa- + &-) 
+n(2apa) — a) tip — pa) a), (2-4) 

incorporates both the driving term £(t) and losses 
through the cavity mirror and by fluorescence. The solu- 
tion p is (/-dependent, and the final result for the density 
matrix is 

p = I P(g)p(g)dg (2.5) 
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wher e p( g) represents the solution of the master equa- 
tion (2.4) for fixed g. 

In order to perform NPCS, which is designed to probe 
the N th couplet, the driving field must be TV-chromatic; 



£{t) = 



N 



m— 1 



(2.6) 



for {£ m } a set of N constants, which are assumed to 
be real without loss of generality. Moreover, the ampli- 
tudes £ m are sufficiently large to ensure significant occu- 
pation of the excited states but not so large that Stark 
shifts or occupation of the higher order (n > N) states 
is significant. By judicious choice of each frequency uj m , 
selective excitation to the N th couplet is possible. 

In the rotating frame w\ = 0, and we define detun- 
ings S m = uj m - uj. 

The randomness of the coupling strength {g} is re- 
sponsible for the inhomogeneous broadening depicted in 
Fig. 0, where g max is the coupling strength between the 
atom and cavity mode at an antinode of the longitudi- 
nal axis. As the probability for the atom-field coupling 
strength being g max is quite small, we choose instead 
to selectively excite the coupled atom-cavity system for 
some other g = gj < <? m ax[j- Thus, we fix lo\ = lu + gj 
as shown in Fig. pL and the amplitude of this field is, of 
course, £\. More conveniently, the normalized detunings 
arc defined by 



<W.9/- 



(2.7) 



We also define a normalised coupling strength as g — 
gj gj- The master equation ([O]) can be expressed instead 
in terms of the Liouville superoperator as 



C(g,t)p(g,t) = Q(g)p(g,t) + £ (£ m e- i ^-^ t E + 

m=2 

-£* e *W»-»/)*E_)p(0,t) (2.8) 

where T,±p = [a±, p], and Q{g) includes conservative and 
dissipative superoperators. The Bloch function method 
is applied by expanding 



(2.9) 



for Z N 1 the set of all length N — 1 vectors with 
integer values and 1 the vector with unity as every 



*The multiphoton spectroscopic signal quality will depend 
on the choices of gf, for our simulation, we choose a value 
which provides a good signal but not necessarily the optimal 
signal. 



component. Transient effects c an be neglected; hence, 
limt-^oo f>k(g,t) —* 0. Thus, eq. (|2.8|) reduces to 



N 

£rn (£ + P£_ fm l (<7) ~ £_P£_ fm+i (s) 



m=2 



tk ■ (S - g f T) + Q(g) p R (g)=0 (2.10) 



where time dependence is ignored as t — > oo, and I m is a 
vector with all elements being except the m th element 
which is one. Writing the superoper ators Q(g) and £± 
as matrices and p(g) as a vector, eq. ( 2. 10] ) represents in- 
finitely many coupled linear equations. In order to reduce 
the number of equations to a finite number, we introduce 
a positive integer q and establish the approximation 



JV-l 

Pk(g) = Vfc satisfying ^ \ki\ > q, 

i=l 

which is valid for sufficiently small {£k}- We perform our 
expansion in the dressed state basis {|n)±|n £ {0} UZ+} 
which satisfies 

H\ri)± =h[m)± Vrig]\n)±, H\0) = 0. 

This spectrum of H eigenvalues is shown in Fig. |l|. For 
NPCS, we truncate beyond the N + 1 couplet. Hence, 
each coefficient pz{g) is of length (1 + 2(N + l)) 2 . Set- 
ting q = 1 leaves 5 + 2(7V — 3) matrix equations, each 
square matrix of dimension (l + 2(7V + l))x(l + 2(iV+l)). 

We are particularly interested in pz—tf{g), which is the 
'dc', or non-oscillating, component of the Bloch expan- 
sion. In an experiment, the iV-chromatic field would 
drive the atom-cavity system, and ./V-quanta resonances 
would yield ./V-photon decays over a timescale shorter 
than the cavity lifetime 113]. To a good approxima- 
tion, the TV-photon coincidence rate is proportional to 
(a^ N a N ), and we evaluate this mean with respect to the 
density matrix component P£_x(g). The elements of the 
density matrix are designated 

p o{g) = (0\Pk = $(g)\0),Pni>(9) = e(n\p R=5 {g)\n') e/ , 
Pon(g) = (o\pk=d(g)\ n )e = PnoG?)*- (2-H) 

The resonance would be observed in practice by fixing Si 
for i — 2, 3, . . . , N — 1, and varying Sn- The iV th cou- 
plet is then observable experimentally as an increase in 
the ./V-photon coincidence rate as a function of Sn', i-e. 
as (a* N a N ) vs Sn. 

However, excited state resonances are complicated by 
the existence of off-resonant excitations which result in 
spurious iV-photon decays. To study these resonant 
and off-resonant effects in greater detail, we calculate 
multiphoton peak heights at various values of Sn using 
the non-Hcrmitian Hamiltonian formalism. The mas- 



ter equation (2.4) incorporates a non-unitary evolution 
which can be treated as a combination of 'loss' terms 
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and 'jump' terms, in the sense of quantum trajectories 
|l0f . The non-Hermitian Hamiltonian is Jl2| 

N 

H eS (g) ^(u-u^Af + igA + iY, £ m (e~ i(tfra_fl/) *0-+ 

m— 1 

_ e i(«m-9/)t ._-) _ Kfl t a _ ( 7/ /2)cr +C r_. (2.12) 

We expand the Hamiltonian as a matrix in the trun- 
cated dressed-state basis, where coefficients of states in 
the (N + l) th couplet and higher are ignored: the time- 
dependent state is approximated by 



\iP(g,t)) « c o (g,t)\0) + c n(9,t)\n) r 



(2.13) 



n—1 s— ± 



with (ip(g,t)\ip(g,t)) < IV t > and co(t = 0) = 
1 and c e n (t = 0) = 0. We work in the rotating- wave 
approximation, and the 2AT+1 coefficients {co, {c^j}} can 
be written as a vector c(g,t). The matrix differential 

equation is c(g,t) = M(g,t)c(g,t). We can write the 
matrix as 



M(g,t) = J2 M te 



■iQ.it 



(2.14) 



where O is a vector of unequally spaced discrete frequen- 
cies, ordered from f2o = to ever-increasing values of 
frequency. That is, Qi > flj for i > j. In solving the 
equation, we ignore terms £le>L where L is a cut-off pa- 
rameter. Physically this corresponds to retaining terms 
responsible for Stark shifts in levels up to order L. Trun- 
cating this expansion is valid because we assume that the 
amplitudes of the driving fields {£ m } are small. 

This calculation of {co, {c^}} allows us to approximate 
the A-photon count rate (NPCR) (a^ N a N ) for given val- 
ues of S 7v where pe aks ar e observed in the full simulation 
following from eq. ( 2.1 0| ) . We use this analysis to verify 
the validity of the computer simulation, applied to the 
special case P(g) = 5{g — <?/), and to observe the impor- 
tance of the Stark effect on the peak heights. As we show 
below for 3PCS, we have excellent agreement between nu- 
merical simulations and this semianalytic approach using 
the non-Hermitian Hamiltonian formalism. 



III. THREE-PHOTON COINCIDENCE 
SPECTROSCOPY 

For large A, significant computer time and memory is 
required to solve the equations, but 3PCS, corresponding 
to N = 3, is readily solved. To probe the third couplet of 
the JC ladder, a trichromatic driving field is employed. 
Ideally a photon of frequency io\ = u> + gf induces the 
transition |0) — > |1) + , followed by another photon of fre- 
quency U02 — oj + (\/2 — l)g/ which takes the excitation 



from the |1)+ to the |2)+ state, and, finally, a third pho- 
ton of frequency 0J3 scans the system over a range of 
frequencies including the |2)+ < — > |3)± transitions as 
shown in Fig. |l|. 



By setting N = 3, cq. (2.1C) reduces to 



0=[i (k 2 8 2 + k 3 8 3 ) + Q(g)} Pk 2 ,k 3 {il) + ^S+Pka-i.kaCs) 
-S 2 'S-Pk 2 +i,k 3 (g) + £3^+Pk 2 ,k 3 -i(9) 
-f 3 S_pfe 2 ,fe3+i(g). (3.1) 

Applying the approximation Pk 2 ,k 3 {d) = for Ifel + lA&l > 
q reduces the number of coupled matrix equations to 
n(q) = 2q 2 + 2q + 1, where coefficients corresponding 
to the fifth couplet and higher are ignored. Hence, there 
are 81 scalar coefficients of Pk 2 ,k 3 (g)- In order to reduce 
computing time, we set q = 1 thereby yielding 2025 si- 
multaneous equations. 

The signature of genuine three-photon decay can be 
obtained by measuring the three-photon count rate 



(3PCR), 



it3A.3\ 



1 S3 to observe three-photon resonance 



peaks. The density matrix elements and the 3PCR are 
shown in Fig. |^ for a range of g. In Fig. |^(a) we observe 
two important features for poo- Firstly, there are two val- 
leys located at S3 = ±g due to the vacuum Rabi splitting 
effect. Secondly, there are two valleys located at g = 1 
and g = v2— 1, independent of S3. The former valley due 
to off-resonant three-photon excitation to the third cou- 
plet from the ground state via the pathway uj± — > uj 2 — > 
L02 (a photon with frequency u>% resonantly excites |0) 
to |1)+, followed by resonant excitation by a photon with 
frequency 101 to |2) + and finally off-resonant excitation 
by a photon with frequency loi to |3)+). The energy of 
the state |3)+ is H(3uj + (2y2— 1)5/) for g = 1. The sum 
of the energies in the three photons producing the ui% — > 
ll>2 —* UJ2 off-resonant excitation pathway is h(3uj+\/3gf ). 
The detuning is thus (2y/2 — 1 — V3)gf O.lgf, which 
is quite small, thus ensuring the significant depletion ob- 
served in Fig. U(a). In Fig.^(b), where p^ is plotted, a 
ridge is observed at g = 1, and the off-resonant excitation 
to yoJ 3 + is thus clear. 

The second valley in Fig. ||(a) occurs at g = \[2 — 1. 
This depletion from the ground state arises due to reso- 
nant excitation from the ground state |0) to the excited 
state |1) + via absorption of a photon of frequency uj 2 
which is fixed (independent of S3). In contradistinction 
to the presence of a ridge at g = 1 in Fig. ||(b), a ridge is 
not observed in Fig. ||(b) at g = \[2 — \ because excitation 
to the second couplet is off- resonant. On the other hand 
excitation to the second couplet is resonant for g = 1 and 
so the ridge is visible for g = 1 in Fig. ||(b). 

The two valleys in Fig. ||(a) at g = 1 and g = \[2 — 
1 have different depths. As each valley is induced by 
resonant excitation from |0) to |1)+, the depletion of poo 
can be calculated from a two-state approximation [pl]Jl4[ 



Poo 



1 



£- 



l ( K + 1I /2r + 2£l 



1,2. 
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For g = 1, m = 1 and £1 = l/y/2 producing p 00 = 13/17 
as observed in Fig. ||(a). Similarly, for g = y/2— 1, m = 2 
and £ 2 — v2 producing p o = 25/41 which again matches 
Fig- 1(a). 

A subtle feature of Fig. g(b) is the presence of two dips 
along the ridge at g = 1. One dip occurs at ^3 = —1. This 
dip is due to competition between two excitation path- 
ways. One path involves resonant excitation via the u>\ 
photon to |1)+, and the second path is resonant exci- 
tation to |1)_ via the U03 photon. Excitation to |1)_ 
diminishes the probability of excitation to |1)+ which is 
necessary for the u>\ — * cj 2 — * LU 2 off-resonant excitation 
to |3)+, hence the dip in p^ . 

The second dip occurs at S3 = — (\/2 + 1) which also 
occurs due to competition between paths; however, this 
dip is less noticeable. Both paths experience resonant 
excitation to |1)+, but the U03 photon excites to |2)_ in 
competition with the u>2 — > ^2 two-photon subsequent 
excitation to |3)+. As the competition occurs for elec- 
trons in the |1)+ state, instead of for the |0) state in the 
case of the other dip, the competition, and hence the dip, 
is less significant. 

Two prominent off-resonant peaks are centred 
about (83, g) = (1,1) and the other at {83, g) = ((y/2 — 
l) 2 , y/2 — 1). The former peak is due to cooperative ex- 
citation pathways via the three-photon excitations u\ — > 
u>2 — * ^2 and u>3 — > 0J2 —* k>2- The equality between UI3 
and wi is responsible for this cooperative effect. The sec- 
ond off-resonant peak is due to cooperation between the 
pathways ll>2 - * ^2 — > ^2 and 0J2 L03 ^ 0J3. The ridge 
at g = \/2 — 1 due to the u>2 — *■ ^2 - ► ^2 pathway is, 
however, negligible due to off-resonant excitation to the 
second couplet as well as to the third couplet. 

There is a prominent peak at (83, g) = (V3- y/2, 1) in 
Fig. ||(a), corresponding to cooperative excitation path- 
ways u>i — > u>2 — > L03 and u>i — > u>3 — » L02 for resonant 
excitation to |3)+. The off-resonant excitation path- 
way u)% — > LO2 — * ^2 also occurs. 

Finally, we observe a peak for g < 0.1 and centred 
at ^3 = 0. This peak is very large near g = (not 
shown) and corresponds to very small splitting of the 
couplets. Consequently, the system behaves much like 
a decoupled atom and cavity and acts as a resonator 
for S3 = 0. The 3PCR, proportional to (aJ 3 a 3 ), is plotted 
in Fig. ||(c). The similarity between Figs. ||(b) and ||(c) 
is evidence that occupation of |3)+ is a good indicator of 
the 3PCR. However, occupation of |3)_ also contributes 
to the 3PCR. The peak at S3 = -(VE+VT) in Fig. §(c) is 
due to the excitation pathway u>i — > 0J2 — > W3 which res- 
onantly excites to |3)_. The peak at S 3 = — (\/3 + \/2) is 
somewhat diminished, however, by the off-resonant u)\ — > 
L02 — > LU2 excitation pathway. The peak at 5 < 0.1 is 
significantly larger in Fig. ||(c) than in Fig. |(b) due to 
contributions from off- resonant excitation to both |3)_ 
and |3)+. 

The peaks, valleys and ridges in Figs. || have been ex- 
plained in terms of excitation pathways. We have intro- 



duced each of these pathways by studying the system in- 
tuitively, but a verification is poss ible using the effective 
non-Hermitian Hamiltonian (2.12) a nd th e approximate 
time-dependent unnormalized state ( |2.13| ) . In Table we 
presented analytical estimates of the 3PCR (a^ 3 a 3 ^ 



est' 



using the non-Hermitian Hamiltonian, and compare to 
the 3PCR, (at 3 a 3 ), observed in Fig. |(c). Although 
the estimates can vary from the observed 3PCR by up 
to a factor of 2.5, the agreement is excellent consider- 
ing that 'jump' terms have been ignored which signifi- 
cantly decrease the peak height. An exception is the dip 
at (£3, (?) = (— 1, 1) where the 'jump' terms are responsi- 
ble for an increase in peak height. Table || provides con- 
firmation of the pathways suggested as being responsible 
for features in Figs. |[ In each case the cut-off parameter 
is L = 1, except for £3 = — (v3 + v2) where L = 2 is due 
to two pathways: one resonantly exciting to |3)_ and the 
other due to off- resonant excitation to |3) + . 

In Figs. ||(a) and |^(b) we present the 3PCR vs S3, av- 
eraged over P(g), assuming a TEMoo mode in the cav- 
ity and a mask for the atomic beam. The expression 
for P(g) is complicated and is provided in the appendix 
of Ref. O. The average 3PCR is thus 



(ot3 3)= Tr(p(g)a^a 3 )P(g)dg. 



(3.2) 



Fig. ||(a) corresponds to a strong coupling of <?//«; = 63. 
This strong coupling ensures that fine features are not 
destroyed by inhomogeneous broadening. By compari- 
son, the case gf/n — 9 is presented in Fig. ||(b), and 
the deleterious effects of inhomogeneous broadening are 
evidently much stronger for the weaker-coupling case. 

In each graph the dotted line corresponds to the 3PCR 
after background subtraction. To perform background 
subtraction, the experiment is performed four times, once 
with all three field on and scanned over S3. The experi- 
ment is then repeated for the three cases (i) £ 1 — £2 = 0, 

(ii) £i^0 and £ 2 = and (iii) £ x = and £ 2 ^ 0. The 
contribution to the 3PCR is then due solely to the ex- 
citations by (i) LJ3 photons, (ii) uj\ and 0J3 photons, and 

(iii) L02 and 103 photons. Case (i) is a subset of (ii) and 
(iii), hence must be subtracted from (ii) and (iii). Using 
the notation of Ref. [jll[, we denote the spectrum after 
subtraction as 

M 3 \5~3) = W3)-FM £l =0(fe) 

-F^ 3 T5 2 =o(^)+F^ 3 T £l = , f2 =o(^). (3-3) 

In Fig. ||, the value of background subtraction is appar- 
ent. The 3PCR for the strong-coupling case gf/n = 63, 
depicted in Fig. ||(a), is improved by background sub- 
traction. Broadening about 83 — is reduced, and the 
peak at 8 = V3 — \[2 is more evident after background 
subtraction. Furthermore, the multiple peak structure 
near 83 = 0, before background subtraction, consists of 
undesirable off-resonant contributions. These are effec- 
tively removed by background subtraction, and the dip 
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in the 3PCR at 63 = 1 is evident. The reason for this dip 
is that LU3 = u>i, and the two pathways u)\ — > L02 — > ll>2 
(which is responsible for the background 3PCR, inde- 
pendent of the value of S3) and L03 — ► L02 — y UJ2 are com- 
plementary. Of course the most important peak occurs 
at 5 = — (y/3 + y/2), which is outside the inhomogeneous 
broadening region. Observing this peak would not re- 
quire the time-consuming background subtraction meth- 
ods necessary for discerning the other peaks. This peak 
at S3 = — (V3 + \/2) exhibits the desired "\/3" signa- 
ture for excitation to the third couplet as well as the 
"\/2" signature arising from the excitation from the sec- 
ond couplet. 

One of the finer features in Fig. ||(a) is the dip at ^3 = 
— 1. This dip is due to the significant dip in Fig. ||(c) 
for (63, g) = (—1,1). However, the size of the dip is 
somewhat reduced due to a very small peak at (63, g) = 
(—1, \f2 — 1) in Fig ||(c) due to the off-resonant excitation 
pathway L02 — > L03 — > lo\ to |3)+. Although this peak is 
quite small, P{g) is more highly weighted for low g. 

For Fig. |^(a) the strong coupling case gf/n = 63 has 
been adopted. The importance of strong coupling is that 
homogeneous broadening due to widths 7j/2 and k are 
small, and contributions due to off-resonant transitions 
are less significant. In previous analyses of PCS, albeit 
for the two- photon case, the coupling strength gf/n — 
9 has been adopted. For 3PCS, such a coupling is too 
small. In Fig. \%b) some of the more dramatic features 
of Fig. ||(a) are still discernible but degraded to a level 
of near indistinguishability from the background. We 
observe the dips at 5 = ±1 and the (very broad) peak 
at 5 — — (vo+ v2)- However, experimental observation 
of such features is unlikely, and a higher coupling strength 
is desirable. A coupling strength of gf/n — 63 is not 
required, but a coupling strength higher than gf/n = 9 
is necessary. 

IV. CONCLUSION 

The technique of photon coincidence spectroscopy 
(PCS) has been extended from driving the atom-cavity 
coupled system by a bichromatic field and measuring 
two-photon coincidences [p|-p^] to multichromatic driv- 
ing fields and multichromatic coincidences. Whereas two- 
photon coincidence spectroscopy (2PCS) enables obser- 
vation of the 2\/2.9 splitting of the spectrum, associated 
with the second couplet, higher-order photon coincidence 
spectroscopy allows direct probing of higher couplets in 
the JC ladder. This technique thus offers a valuable tool 
for resolving higher-order spectral phenomena in cavity 
quantum electrodynamics. Although the focus here has 
been on an atomic beam and the spectrum associated 
with the JC ladder, this scheme can be adapted to study- 
ing spectra of other cavity QED systems with a discrete 
spectrum and significant inhomogeneous broadening. 



We can also see that off-resonant phenomena are in- 
creasingly important for probing higher-order couplets. 
The background subtraction scheme is more intricate. 
Also resolving peaks for higher-order couplets requires 
increasingly large coupling strengths. The required cou- 
pling strength places a bound on the feasibility of TV- 
photon coincidence spectroscopy, and we see that a large 
coupling strength is required even for N = 3. 
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FIG. 1. A three-photon excitation scheme from the 
ground state |0) to the first three couplets \n) e (n < 3) of 
the dressed states with inhomogeneously broadened energy 
bands. 

FIG. 2. Plots of (a) poo, (b) p\t and (c) the (a t3 a 3 ) 
vs normalised scanning field frequency, S3, and coupling 
strength, g, for Ji/k = 1, = 63, £i/k = 1/V% and 

E'ljn — £-j,/k — \/2. 



FIG. 3. Plots of (a t3 a 3 ) with (solid) and without (dots) 
background subtraction vs 83 /it for Ji/k = 1, £i/k = l/y/2, 
£2/1* = £%Ik — \/2 averaged over P(g) for randomly-placed 
atom in the TEMoo mode for (a) g; j n = 63 and (b) g/ /k = 9. 
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S3 


(at ; V) 




(V2 - 1) 




1.2 x 1(T 3 


3.1 x 1(T 3 


1 


1 


1.7 x 1(T 3 


1.7 x 1(T 3 


1 


V3-V2 


1.5 x 1(T 3 


2.6 x 1(T 3 


1 


-1 


2.4 x ifr 4 


2.1 x 1(T 4 


1 


-(V2+1) 


3.3 x irr 4 


3.3 x 1(T 4 


1 


-(V3 + V2) 


1.3 x 10~ 3 


2.1 x 1(T 3 



TABLE I. Three-photon count rate (3PCR) (a t3 a 3 ) 
and approximate 3PCR (a t3 a 3 ) ost calculated via the 
non-Hermitian Hamiltonian formalism. 
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